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I. INTRODUCTION 


In the past, passive sonar systems were only able to provide bearing information, and 
the estimation of source range relied on analytical, over-simplified techniques. Starting in 
the mid-60’s, however, with relatively advanced computational resources, the received 
signals have been processed in order to also provide information about the distance to the 
source. Today, one of the most used and successful research techniques is matched-field 
processing (MFP) where a comparison is made between the parameters of the received 
signal and those of signals generated by a synthetic source virtually positioned at each point 
in the search grid. A match is obtained when the correlation of the synthetic and the true 
signals are at a maximum. 

Transient signals are especially important for naval purposes, since they can represent 
several instances of a ship’s routine, such as firing a torpedo, opening of hatches, starting 
of pumps, etc. Therefore, the Transient Localization Project at the Naval Postgraduate 
School has been studying, since 1993, several algorithms for localizing transient sources 
based on methods similar to MFP. 

In 1995, three distinct schemes were studied: (a) the fully coherent localization, which 
was considered impractical as it required an extremely precise (and unavailable) 
propagation model, along with accurate environmental descriptions; (b) the semi-coherent 
localization, which correlated only the slower varying amplitude terms of the signal, 


trying to surpass the limitations of method (a); and (c) the autocorrelation matching, which 


compared the autocorrelations of the received and modeled signals using a simple Bartlett 
technique. The NPS results showed that method (b) was not successful in localizing the 
signals of naval interest, and method (c), while not perfectly consistent, was the most 
useful technique. According to Miller et al. (1996), autocorrelation matching has an 
advantage over the other methods, since its processing allows the notching of the noise 


energy and the possibility of a multi-transient gain. 


A. THE NUWC-TSP EXPERIMENT 

The transient signal processing (TSP) experiment was a field study ined by 
NUWC, and conducted both by NUWC and C&M Tech., Inc. on 28-31 July 1997. Its 
purpose was to provide data to study the influence of a highly variable shallow water 
environment on acoustic propagation, and the effectiveness of localization algorithms. 

The study area is shown in Figure 1.1 and as it is influenced by oceanographic and 
geological processes related to the shelfbreak, it provides a complex coastal environment 
for sound propagation. The continental shelf water is cold, and it presents large variations 
with location and season. The slope water, on the other hand, is warmer and influenced by 


the Gulf Stream.(Miller, C. W., 1998) 
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Figure 1.1. The NUWC-TSP Experiment Geographic Location. 


Synthetic signals were generated and transmitted in a variety of scenarios which were 
differentiated by the source depth, source-receiver range, and bearing. Two projectors 
(ITC-104 and HX: 188), a 32-element vertical-line hydrophone array (VLA), and data- 
recording equipment constituted the main assets of the experiment. The VLA 
configuration is depicted in Figure 1.2, and the experiment configuration for one of the 


scenarios is shown in Figure 1.3. 
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Figure 1.2. The VLA Configuration. 





Figure 1.3. The NUWC-TSP Experiment Configuration. 


It is important to note that hydrophones 17 and 18 were non-operative throughout the 
experiment. In this work, we were only interested in the linear frequency modulated 
(LFM), short duration (100 ms) signals - the transmitted chirp and a sample of a received 
burst can be seen in Figure 1.4. Their frequency spectra are observed in Figure 1.5. 

The processed data corresponded to three particular cases: scenario 9 (source located 
at 4650 m range, 45.7 m depth), scenario 11 (source at 460 m range, 45.7 m depth), and 
scenario 16 (source at 2194 m range, 27.4 m depth). The information on source location 
was provided by NUWC, based on ship’s log and GPS data. While some uncertainties are 
present in this information, inherited from the processes used to gather it, we assumed it as 
our reference. The recorded tapes also contain array tilt information, which was neglected 


in this work for simplicity. 





Figure 1.4. The Transmitted LFM Burst (Left) and the Received Signal at VLA (Right) in 
the Time-Domain. Note That the Axis Are Different on Each Plot. 
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Figure 1.5. Frequency Spectra of Transmitted Chirp (Left) and Received Burst (Right). 


In order to support the analysis of the experimental data, six conductivity-temperature- 
depth (CTD) casts were made. Two of them did not present useful data and were 
disregarded. The other four provided the sound speed profiles (SSP) shown in Figure 1.6, 
which also displays the date and time (Zulu) they were taken. For processing, however, we 
used only SVj97_08 (applied to scenario 9), and SVj97_12 (applied to scenarios 11 and 
16), because they were gathered at times closer to those of the studied runs. Still, it is 
important to note the large variations, particularly near the lower part of the thermocline 
‘between 20-35 m. The associated pycnocline supports the propagation of nonlinear soliton 
waves which travel from the edge of the slope onto the shelf. Such oceanographic features 
can create localized regions of high sound speed contrast as strong as 25 m/s over a 100 m 


horizontal range in the direction of the soliton propagation (Chiu et al., 1997). 
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Figure 1.6. Sound Speed Profiles (SSP’s) Collected During Experment. 


B. THESIS OBJ ECTIVE APPROACH AND OUTLINE 

The objective of this thesis is to examine the influence of shallow water acoustic 
variability on frequency-domain autocorrelation matching (FACM) localization. In order 
to perform such analysis, we processed the acoustic VLA data, identified the LFM bursts, 
and high-pass filtered them to reduce the array strum influence (cut-off at 50 Hz). The 
resulting bursts, as well as the replicas generated by the MMPE propagation model, were 


applied to a Bartlett-type processor based on matching broadband signal autocorrelation 


localization algorithm). The acoustics variability study followed with the plane-wave 
beamforming, modal amplitude structure, and spatial coherence analyses. 

The remainder of this thesis consists of five chapters. Chapter II describes the 
Monterey-Miami Parabolic Equation (MMPE) propagation model. Chapter II discusses 
the Frequency-domain Autocorrelation Matching (FACM) localization algorithm. Chapter 
IV describes the techniques used to evaluate the environmental Sfkences and process the 
signals. Chapter V presents the localization results. Chapter VI concludes this study with a 


summary of findings and recommendations for future work. 


Il. THE MONTEREY-MIAMI PARABOLIC EQUATION 
PROPAGATION MODEL 


The Monterey-Miami Parabolic equation (MMPE) acoustic propagation model 
(Smith, K., 1996) was used to generate the replicas which simulated the received signals 
from a synthetic source on the gridpoints of a 2-D search space (depth x range). It allowed 
us to predict the arrival structure at the VLA due to a transient-like, broadband, point 
source. In this chapter, we introduce the general theory behind the parabolic equation 
model, as well as the method used for its implementation - the Split-Step Fourier (SSF) 
method (Hardin and Tappert, 1973). 

The parabolic equation method for computing underwater acoustic propagation was 
introduced by Tappert (1977). To begin the parabolic approximation development, 


consider the representation of the time-harmonic acoustic field in cylindrical coordinates 


-i2nft 


P(r, 2,0,f,t) = pAr,z O)e (1) 


As the model is based on an approximation to the Helmholtz wave equation, substituting 


Eq. 1 into the latter in cylindrical coordinates leads to (Jensen et al., 1994) 


2 
12/5 19D, O Pr 


"5 r= 2092 * > a +k2n? Pr = _4nP JO(r—rs) . (2) 


a 


c(7, z, 0) 


where n(r, z,0)= is the acoustic index of refraction, c, is the reference sound 


O. . 
speed, c(r, z, O) is the acoustic sound speed, and k,= — is the reference wave number. 
C 
O 


In this derivation, however, we neglected density variations, which could be 
incorporated into a new index of refraction without any loss of generality. Note also that 
the environmental characteristics are within c(r,z,), and the reference pressure level 


P is defined as the pressure amplitude at a distance R, = 1. 


By assuming the ocean acts as a waveguide with a cylindrical coordinate system, 
acoustic energy 1s primarily propagated outward from a source in the horizontal direction. 


The pressure field, therefore, can be approximated by 
PAr. zo) = Wr, z, OHS) (kyr), (3) 
where the slowly varying envelope wAr,z,) modulates the outgoing zero-th order 


Hankel function of first kind H{)(k,r). 


In the far field (A,r » 1), we can approximate the Hankel function’by the first term of 


its asymptotic expansion (Gradshteyn and Ryzhik, 1994), 


, 2 Boe (4) 
imk r 
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AW (k r) = 





and Eq. 3 becomes 


R 
PAT, 2,0) = p. fey (7.2.0) (5) 
Note that Eq. 5 is normalized such that at r = R,, |w| = 1 and [pq = P,. It expresses 


the relationship between the PE field function w(r,z,) and the acoustic pressure 
pdr, z, >). By substituting it into Eq. 2 and dropping the source term on the night hand 
side, we obtain 

2 2 


If we neglect the azimuthal coupling and the far field terms, Eq. 6 can be written as 


Loy _ 
kor ae Cree onc (7) 
where 
_ isan 
Top ™ 55. (555) ® 
and 
eb (9) 
Oop” 2 ; 
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Note that this is only valid for re «WwW, once we used the uncoupled azimuth 


approximation. 


Eq. 7 is known as the standard parabolic equation (SPE), and the accurate solutions are 
limited to a half beam width of 10° to 20° for the propagation angle. However, in order to 
extend this limit to 40° ~ 70° , a higher order wide-angle parabolic equation (WAPE) 
approximation (Thomson and Chapman, 1983) can be used. It also features less sensitivity 
to the choice of k,, and to the phase errors in typical deep ocean conditions, relative to the 


SPE approximation.(Jensen et al., 1994 and Chin-Bing et al., 1993) 


The WAPE operators are defined as 


a 
a. | ae 
[eee alt + ie is 1 (10) 
and 
Ue =U. ieee (11) 


In order to numerically solve the parabolic equation, the MMPE uses the split-step Fourier 


(SSF) method. This algorithm integrates the solution in range by applying the operator 


Uyaspr in the z-domain., and the operator Ty, pz in the k, -domain. The latter is defined, 


in the k, -domain, as 


~ k2 
TwaPE(k,) =1- fn -#). (12) 
Oo 


Note that both operators are just scalar multipliers, and may be applied independently. 
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The approximate solution for yw is iterated in range according to the following 


expression 


~ik,ArUwape(r, Z) —ik,ArTwape(r, kz) 


w(r+Ar,z) = e FFT ¢ [FPT(W*(r,2))I*), GES) 


where the FFT calculation is consistent with the form described in (Press et al., 1988). A 
forward FFT and double conjugation emulate the inverse FFT. The model outputs the field 
functions y (both magnitude and phase), referenced to a unit magnitude at r = Im, and 
computed at the spatial grid points. 

The source assumed in our model had a 600 Hz bandwidth centered at 900 Hz. 
Therefore, the acoustic field was computed for each one of the discrete frequencies 
(N=512) of the considered BW, and the model properly represented the broadband 
acoustic field. | 

As the single component of the general field can be expressed as Eq. 5, and assuming 
a windowed (Hanning), normalized source amplitude, the time-domain complex pressure 
field can be written as | 


oo 


—i2nft 


p(r,z.t)= |] plr.zfle af 
il a -i2nft 
= =) ee Ne df . (14) 
Lr —oo 


i2nf— 
vice Cc. : : 
Note that the overall phase factore ~ = e can be neglected if we consider that the 


arrival times are given in terms of “reduced time” 7 = t— | , then 
O 


pul z) ve, z, fre af. Gls) 


e ° ° ° a 
This means that the time domain is heterodyned around t = —, and the use of a reduced 
O 


time does not have any influence on the autocorrelation function. The model also 
heterodynes the signal by shifting the center frequency to zero (d.c.), in order to reduce the 
computational load associated with large transform sizes. This procedure does not 


introduce any consequence to our algorithm. 


The MMPE model, a FORTRAN® code, is the latest version of what was formerly 
known as the University of Miami Parabolic Equation (UMPE) model (Smith and Tappert, 
1993). It requires a series of input files with strict formats, which properly represent: 

¢ environmental data (sound-speed profile, bathymetry of water/bottom interface, 

bathymetry of the deep layer beneath the water/sediment interface, acoustic 
parameters of this “deep” layer, and acoustic parameters of the medium just below 
the water/bottom interface); 

¢ source data, with two types available, wide-angle and vertical line array (steering 

allowed). We only used the wide-angle, which approximates the point source. The 
center frequency, bandwidth, and number of discrete frequencies (a power of two, 
because of the FFT) are also required. 


All these input files are specified in a main one, the “pefiles.inp’, where several other data 


(name of output files, grid size and grid steps, Cg, size of vertical FFT, etc.) are contained. 


The output is presented in a single binary file composed of a header (with the information 
needed for post-processing), and the complex PE field function y at selected grid points 


for every discrete frequency. 





Ill. THE FREQUENCY-DOMAIN AUTOCORRECTION MATCHING 
(FACM) LOCALIZATION ALGORITHM 


The algorithms aia in the Transient Localization Project at Naval Postgraduate 
School are designed for data from a single receiver hydrophone and a point source. The 
basic concepts of the correlation functions are outlined in Bendat and Piersol (1971), and 
the routines are described in Miller et al. (1996), Pierce (1996), and Hager (1997). From 


Brune (1998) we have that: 


‘Localization algorithms may be considered generalized 
beamformers in which the plane wave replicas have been replaced by 
more complicated replicas of the acoustic propagation (e.g., modes, 
beams, or the vertical pressure field). The algorithms, usually referred to 
as processors, are in most cases based on a Hermitian quadratic product. 
The exact form is determined by the constraints that are put on the 
processor output. . 


As discussed in Chapter II, replicas are simulated received signals of synthetic sources 
located at the points of a search grid. Generated by a propagation model, they are matched 
with the received signal. Considering that the source functions are equal and a perfect 
model is used, the exact match should occur where real and synthetic source locations 
coincide. An ambiguity surface usually represents the localization algorithm results. 

The reciprocity principle states that, in a time-invariant and homogeneous 
environment, the acoustic pressure at location A due to an omnidirectional source located 
in B is equivalent to that of location B due to an omnidirectional source located in A. 


Using this principle, we reduce the computational load required to obtain localization, 


since we are able to set the propagation model to generate replicas from the receiver 
position to all source possible locations in the grid. 

In this chapter, we discuss the derivation of the frequency-domain autocorrelation 
matching (FACM) localization algorithm. Autocorrelation matching was originally 
designed for use in the time-domain (TACM), therefore it is more convenient to start its 
study treating the problem in time, and then to convert the results to the frequency domain. 
Previous investigations done by de Kooter (1997) showed that the TACM algorithms 
present very small footprints and large phase errors at higher frequencies, being useful for 
very low frequency signals. It was also observed that matching the time-domain 
amplitude-squared signals, which seemed relatively similar, could be more effective and 
robust. Therefore, as we will describe in the following paragraphs, an almost natural 
further development led to the frequency-domain algorithm. We begin with a complex 
pressure time series P(t), which represents a detected transient arrival. 


The signal autocorrelation in the time-domain can be expressed by 


Tpp(t) = | P (t)P(t+t)at (16) 
or, in terms of the frequency-domain response, by 


—i20fTt 


Tpp(t) =| P (Phe af, (17) 


where 


P(f) = fu P(e dt. (18) 


Then, the normalized autocorrelation becomes 


Tpp(t 
Cpp(t) = — pe | ; (19) 
PON af 
A similar derivation apphed to the predicted replica signal R(t) leads to 
Tpp(T) 
Cpa(t) = (20) 
{ lrol'a 


The autocorrelation matching algorithm is based upon the inner product of these two 


quantities and, normalized, in time-domain, it can be expressed by 


| GG ee ac 
a a a TT 21) 
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By analogy, in the frequency-domain it becomes 
; J CeR(OC5p(O) a0 
ARP = So (22) 
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where C;, and C;, are the autocorrelations in frequency-domain of the transient signal 


and of the replica, respectively. 


Since 
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Eq. 22 can be rewritten as 
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In order to reduce the influence of noise, we shall remove the zero-lag component of 
the autocorrelation function, which is equivalent to removing the mean of the squared 


amplitude of the time-domain signal. Therefore, we define 


IR'(d)|? = IRD - (IRC): (25) 
and 


, 2. 2 2 
IPA = IPCI — CIPO): (26) 
where < >, is an average over all times. Eq. 24 then becomes 


{RPP oat 
ARE a ee Sey ey 
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which is focused on the influence of matching the non-zero lag components. 


~~ 
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However, it still presents an ambiguity in (absolute) tme which was not present in 
TACM, and we shall introduce a search parameter T to allow us to slide the replica signal 


in time, searching for the optimal match. The FACM function then becomes 


J IR+ oP rar 
ARP = max, = 


SS 5 arr (28) 
| IR’(t)| ‘de Preoltar] 


where max, is the maximum value of the function for any value of the search parameter 


A MATLAB ® implementation of the algorithm based on Eq. 28 was used to generate 


all the ambiguity surfaces shown in this work. The sliding-t operation was performed 


according to the following development. Let g(t) =|P'(t)|’ and h(t) =|R'(1)|", then 


si =f andr, (29) 
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and 
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Therefore, the numerator of Eq. 28 becomes 


oO 


Alec) = [ n(e+ c)g(a)ae. | | (31) 


If A(t) 1s a real function, 


oO 


A(T) = [ n*(0+ De(nat 

= | [ PH (He af [eG af’ dt 
= iL H* (Ge | peer — St df’ df 
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= | H*NGHe df , (32) 
knowing that 

fe "ae = 86° =f). (33) 
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Hence, the search over T was in fact computed by a multiplication in the frequency 


domain and an inverse FFI operation represented by the numerator A(t). Prior 


investigations (Brune, 1998 and de Kooter, 1997) had shown that the FACM algorithm is 


better suited for larger bandwidths and higher frequencies, when compared to its analog in 


ae 


the time domain, the TACM algorithm. Also, FACM is more robust against environmental 


mismatch and presents better peak-to-sidelobe levels and larger footprints. 
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IV. SIGNAL PROCESSING AND ENVIRONMENTAL 
ASSESSMENTS 


A. MATCHED-FILTER PROCESSING 

Matched filtering (MF) is a process where the received signal r(t) is correlated against 
a replica of the transmitted signal s(t), in order to obtain a better signal-to-noise ratio 
(SNR) and reduce the effective temporal structure of the source down to an ideal, coherent 
pulse. It is also known.as a conjugate filter, since its frequency response 1S basically the 
conjugate of the (transmitted) signal spectrum (Tolstoy, A., 1993). In this work, matched- 
filtering was used first to identify the short duration bursts 0.15) in the 2-min. long data 
files. Second, but more importantly, to process the bursts and obtain data emulating the 
existence of a coherent source, instead of the actual non-coherent one. This also provides a 
baseline for the localization results since the model-generated replicas are based on a 
coherent source. Therefore, this is equivalent to localizing a transient with known source 


function. 


The process was implemented digitally in MATLAB ® , according to the expression 


i2nft 


pit) =f S*(NRDe? af. 34) 


where p(t) is the matched-filtered output in time-domain, S*(f) is the complex conjugate 


of the transmitted signal in the frequency domain, and R(/) is the received signal in the 


frequency domain. 
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An example of the MF output in time domain can be visualized in Figure 4.1, and its 
frequency spectrum in Figure 4.2. The raw data is also shown and, by comparison, we can 


observe the significant noise reduction. 


Maichedtite ed burl - Hyd @ 22, Scenano 9 





Figure 4.1. Examples of a Raw (Left) and an MF Burst (Right) in the Time-Domain. 
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Figure 4.2. The Frequency Spectra of a Raw Burst (Left) and of an MF Burst (Right). 
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B. PLANE-WAVE BEAMFORMING 

Plane-wave beamforming is one of the many techniques used to distinguish particular 
features of a received signal while reducing the influence of noise - which is usually 
generated near the surface, and tends to occur in high modes or high angles of propagation 
(Jensen et al., 1994). It consists of decomposing the signal into plane waves by summing 
the outputs of spatially distributed sensors - in this work, eight hydrophones for each of 
the three sub-arrays in the VLA. 

In underwater acoustic propagation, much of the information about the environment 
and the relative source/receiver location is contained in the vertical structure of the 
acoustic field. Therefore, the ability to resolve the vertical structure of the arrival at the 
receiver can lead to estimations of the corresponding source location. 

A beamformer is especially useful because it allows us to give preference to one 
direction of propagation over another, implementing a spatial filter. In this work, it was 
used to observe the dominant beams when the signal was arriving at the VLA, providing 
information for a comparison of the localization algorithm results by the distinct arrival 
angle structures. 

An excellent discussion about the beamforming process is presented by Defatta et al., 
(1988) and its summary follows. Consider an infinite, homogeneous medium containing a 
line array of equally spaced elements positioned along the y-axis, and an undetermined 


number of remote sources. The output of an element located at the origin of coordinates 


Zt 


due to the / source is defined as S )(t). Assuming plane-wave propagation in the far-field, 


a source from direction 0, outputs 


e,(t) = (t+ Gy 


ndsin(8,) 
ara 


where nis the array element index, d is the spacing of the elements, 9 j is the angle of 


arrival from the J” source, and c is the speed of propagation. To accomplish beamforming, 


we must apply weights 0 and time delays to the individual element signals, and then 


coherently sum them. Therefore, the output of the beamformer in the direction 2 1S 


N- 
n= Q 


nd{sin(@,) ‘ea 


(36) 
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To derive the directional response characteristics of the beamformer, we ‘shall consider a 


continuous, complex plane-wave signal propagating across the array, 


iw,t 


Scan | Gi) 


Due to the arrival angle 0 I? the sensor outputs 


1@),t ink,dsin(@;) 


e,(t) =e e ' (38) 


In «= 


where ky = 1 = — js the wave number of the I” source. 
I Es 
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The beamformer output is then 


i= 
b (t)= Pg c TC sa w(8,)e (39) 
a) 
where 
N- 1 
we) = ¥ wy eiidlsin( 61) ~ sin(@n) (40) 
n= 0 


is the spatial Fourier transform of the array weights. 


The beamforming operation applied in this work was implemented using FFT in 


MATLAB®. Single frequencies considered, the time-delay operations described above 
are equivalent to phase shifts that are linear functions of the element index. FFT 
beamforming is the implementation of these phase shifts through FFT operations. The 
processed results for all three sub-arrays, and the three distinct scenarios can be seen in 
Figures 4.3, 4.4, and 4.5. They describe the physical angle of arrival as a function of the 
reduced time, and provide an easy visualization of the multipath nature of the burst 
propagation. The results for sub-array 2, however, are noticeably affected by the lack of 
data from hydrophones 17 and 18. We can verify in Figure 4.3 (scenario 11, range ~500 
m) that the energy is not as concentrated in the lower angles as in the other two figures 
corresponding to larger ranges. However, in all of them we observe the early arrivals of 
the lower angles, the short period in which the energy is contained, and the distinct 


multipath structure. 


29 


Tranerrission Loss (dB re 1m) ~- Scenario 11, Sub—array 3 








Amval Ange (deg) 


0.46 0.46 0.6 052 0.54 0.58 0.58 
Time (ee@0)} 


(a) Sub-array 3 


Tranemesion Loss (dB re 1m) - Scenario 11, Sub—avray 2 


| 





Artal Angte (deg) 


0.46 0.48 0.6 0.54 0.66 0.66 


(b) Sub-array 2 


Tranemesion Loes (dB re 1m) — Scenario 11, Sub—asray tf 


Amval Angle (deg) 





(c) Sub-array | 


Figure 4.3. Beamformer Outputs for Scenario 11 (Range 460 m, Depth 45.7 m). 
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(b) Sub-array 2 
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Figure 4.4. Beamformer Outputs for Scenariol6 (Range 2194 m, Depth 27.4 m). 


3] 





Transmission Lose (dB re im) — Scenario 9, Sub—array 3 


Amval Ange (deg) 





0.52 0.54 0.66 0.66 0.6 0.62 0.64 0.66 0.66 
Time (sac) 


(a) Sub-array 3 
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(b) Sub-array 2 
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Figure 4.5. Beamformer Outputs for Scenario 9 (Range 4650 m, Depth 45.7 m). 
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C. MODAL ANALYSIS 

Normal modes are another physical, frequency-dependent representation of the 
underwater acoustic wavefield. For a given sound channel, the acoustic modes represent a 
unique, natural set of standing waves oscillating along the depth axis (Chiu and Ehret, 
1994). 

The method involves solving a depth-dependent equation, and the complete acoustic 
field is then constructed by summing up contributions of each of the modes weighted in 


accordance to the source depth.The modal equation 





2 
“| 1 au) ® 2 
a th —k rr oe 4] 
pC x Ou ey ™ me (41) 
with boundary conditions ‘Y(0)= O (pressure-release surface at z=0) , and a =) 
z=D 


(perfectly rigid bottom at z = D) is a classical Sturm-Liouville eigenvalue problem, whose 


properties are well known (if we assume that p(z) and c(z) are real). Some are: 


¢ the modal equation has an infinite number of solutions which are like the modes of 
a vibrating string; 


¢ the modes are characterized by a mode shape function ‘YY, (z) (an 
eigenfunction) and a horizontal propagation constant krm (eigenvalues). These 
constants are all distinct and analogous to a frequency of vibration; and 


¢ the modes are unique (orthogonal) and form a complete set, which means that we 
can represent an arbitrary function as a sum of the normal modes. 


BD 


Based upon these and through a modal decomposition implemented in MATLAB, we 
were able to quantify the amplitudes of the first 30 modes as functions of the signal 
(discrete) frequencies and, by IFFT, of the reduced time. The SSP’s used were the ones 
discussed in Chapter I and depicted in Figure 1.6, applied to their respective cases. A 
visualization of the first ten modes of the assumed waveguide for a frequency f = 900 Hz 
is presented in Figure 4.6. The final results of the modal decomposition of the received 
signal for each studied scenario are shown in Figures 4.7, 4.8 and 4.9. Note the energy 
concentration in the lower modes, especially scenarios 9 (5000 m) and 16 (2500 m). This 
is consistent with stripping of higher modes at longer ranges. The energy observed in the 
higher modes must be treated with skepticism since the array does not sample those modes 
very well. The lack of tilt information in this analysis also makes the higher mode 


decomposition suspect. 


Mode #1 Mode #2 Mode #3 Mode #4 Mode #5 Mode #6 Mode #7 Mode #8 Mode #9 Mode #10 
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Figure 4.6. The First 10 Normal Modes (f = 900 Hz). The Bottom Depth is 89.9 m. 
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Figure 4.7. Modal Amplitudes vs. Reduced Time - Scen.11 (Range 460 m, Depth 45.7 m). 
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Figure 4.8. Modal Amplitudes vs. Reduced Time-Scen. 16 (Range 2194 m, Depth 27.4 m). 
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Figure 4.9. Modal Amplitudes vs. Reduced Time-Scen.9 (Range 4650 m, Depth 45.7 m). 


36 





D. SPATIAL COHERENCE (CORRELATION) 

In order to analyze the spatial coherence of the signal’s vertical structure at the VLA, 
single records from each of the three scenarios were considered. We computed the peak 
correlation of each phone’s data with a reference phone (hyd #1, closer to the bottom; 16, 
mid-array; and 32, closer to the surface) normalized such that the autocorrelation of the 
reference phone data was unity. The results provided eee about how unique the 
signal was across the space (vertical direction). We can observe in Figure 4.12, from 
scenario 9 (range ~5000 m), that the signal correlation decreases to values around 0.5 over 
about 5 meters in depth for the reference hydrophones 1 (bottom) and 32 (top), which 
corresponds to one element spacing and approximately two to four wavelengths. The same 
effect can be verified in Figure 4.11 which displays data from scenario 16 (range ~2500 
ey) These results suggest that distinct multipath arrivals at any given time are interacting 
at different locations over the array length. For scenario 11 (range ~500 m) runs, however, 
the cross-correlation is not as low for the same difference in depth, as shown in Figure 
4.10. This is a consequence of the short distance between source and receiver which does 
not allow large variations in the signal nor many multipath arrivals. Also, by observing 
Figure 4.12, we notice that hydrophone 16’s signals are more correlated than the ones 
from hydrophones 1 and 32 for scenario 9. This could be a consequence of its relative 


position (mid-array), but the same effect was not verified in scenarios 11 and 16. 
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Figure 4.10. Phone-to-Phone Peak Correlation - Scenario 11 (Range 460 m). 


38 


Phane-to—phone peak coneiation for Sc enano 16 — MF 


Peak correlation (normalzed) re Hyd 32 
2 2 2 ° 2 2 ° 
w& a wn a ~ © © =. 


o 
bs 


8 


3S 40 45 65 79 78 a0 


SS 
Depth [rm] 


(a) Reference Hydrophone # 32 (Depth 28 m) 


Phone-to—phone peak corretation for Scenano 16 — MF 


Peak caraation (normalized) re Hyd 16 
9° 
uw 


30 3s 40 4S 6S 70 7S 80 


SS 
Depth {m) 


(b) Reference Hydrophone # 16 (Depth 54.8 m) 


PRone-to—phane peak correlation for Sconano 16 - MF 
1 


o9 


Peak conmealion (normalized) re Hyd 1 
° =) ° ° ° 2 ° ° 
+ by bo a uw r.) ~N © 


° 


30 3S 40 4s SO SS 60 6S 70 75 a0 
Depah [mm] 


(c) Reference Hydrophone # 1 (Depth 82 m) 


Figure 4.11. Phone-to-Phone Peak Correlation - Scenario 16 (Range 2194 m). 
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Figure 4.12. Phone-to-Phone Peak Correlation - Scenario 9 (Range 4960 m). 


V. LOCALIZATION RESULTS 


This chapter presents a compilation of results obtained after a series of runs of the 
FACM algorithm. For every scenario, distinct receiver depths were tested by selecting the 
hydrophone data to be used. 

Regarding the environmental parameters used as inputs to the MMPE model, only the 
sound speed profile varied according to the respective scenario. We used only the 
bathymetry along the North radial, which is an approximate 89.9 m isobath with bottom 
composed mainly of medium to coarse grained sands. This type of sand has the following 
properties: compressional sound speed of ~1800 m/s, sound speed gradient of 25s"!, shear 
speed of ~250 m/s, density ~2.03 g/em?, and compressional and shear wave attenuation 
coefficients about 0.1 dB/m and 1.0 dB/m, respectively (Smith et al., 1998). During the 
experiment, the data was recorded at a sampling frequency of 4934.0 Hz. For processing, 
we extracted the received hydrophone response to the 0.1s bursts in files of 4210 points 
which corresponded to 0.8533 seconds. The reason for this procedure is that, in the 


MMPE model output, we have 512 discrete frequencies over a 600 Hz bandwidth. 
Therefore, Af = 1.171875 Hz, and we match T = re = 0.8533. After, we processed it 


again to reduce its length to 2048 spectral points over a bandwidth of 600 Hz to 1200 Hz. 


The model output was zero-padded (both ends) to also contain 2048 points. 


4} 


The spatial grid was, in most of the runs, 128 points in depth vs.:300 points in range. In 
general, it corresponded to a resolution of 0.78 m in depth and 10 m in range. For all cases, 
the source/receiver relative position proved to be a very important factor to the quality of 
the results. This can be observed in Figures 5.1 (a) and (b), where the ambiguity surfaces 
for scenario 9, based on hydrophone 16 (a), and hydrophone 2 (b), are shown. The source 
is located at depth 45.7 m, range 4650 m, and its location is depicted as a white dot. While 
the algorithm with hydrophone 16 (depth = 54.8 m) data was able to localize it within the 
some tolerance, the run with hydrophone 2 (depth = 76.3 m) data did not provide 
localization at all. 

By processing the matched filtered signal, we simulate the (idealistic) situation of a 
coherent source, and it improved the localization results a little. However, even when 
processing the raw data, we obtained satisfactory results where the main differences were 
the bigger footprint and the smaller values of the ambiguity surface. These differences can 
be observed in Figures 5.2, 5.3, and 5.4, which display results for scenario 9, scenario 16, 
and scenario t1, respectively. 

Observing Figures 5.3(b) and 5.4(b), where the ambiguity surfaces from runs with MEF 
data from scenarios 16 and 11 are shown, respectively, we verify that the FACM 
localization algorithm performed reasonably well (errors < 11% in range, and < 40% in 
depth). It is important to note that the sound speed profiles used were gathered four hours 


apart from experiments of scenario 9, more than two hours apart from scenario 11’s, and 
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more than 10 hours apart from scenario 16’s. The mismatch between the actual and the 
modeled environmental parameters contained in the SSP is likely significant. 

We also observed that the localization error decreased with increasing range. This was 
expected since the autocorrelation techniques are strongly dependent on the uniqueness of 
the multipath structure. The error obtained for scenario 9 (range ~5000 m) runs were in the 
order of 2% in range and 30% in depth. However, scenario 11 (range ~500 m) runs 


presented errors as large as 20% in range and 40% in depth. 
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Figure 5.1. Ambiguity Surfaces for Scenario 9-(a) with Raw Data from Hyd.16 (Depth 
54.8 m), and (b) with Raw Data from Hyd.2 (Depth 76.3 m). The White Dot Depicts the 


Actual Source Location. 
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Figure 5.2. Ambiguity Surfaces for Scenario 9 - (a) with Raw Data and (b) with MF Data 
from Hyd.16 (Depth 54.8 m).The White Dot Depicts the Actual Source Location. 
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Ambiguity Surface - Raw data, Hyd.#16, Scenario 16 
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Figure 5.3. Ambiguity Surfaces for Scenario 16- (a) with Raw Data and (b) with MF Data, 
both from Hyd. 16 (Depth 54.8 m).The White Dot Depicts the Actual Source Location. 
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Ambiguity Surface - Raw data, Hyd.#16, Scenario 11 
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Figure 5.4. Ambiguity Surfaces for Scenario 11- (a) with Raw Data and (b) with MF Data, 
both from Hyd. 16 (Depth 54.8 m).The White Dot Depicts the Actual Source Location. 
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VI. CONCLUSIONS 


A. SUMMARY OF FINDINGS 

Several techniques helped us to assess the properties of the acoustic field in which we 
were applying the localization algorithm and infer their influences on it. Through spatial 
coherence analysis, we were able to verify that more unique signals led to better 
localization results, therefore, ratifying Brune (1998). Analyzing the beamformer outputs, 
we can affirm that more stratified, discrete structures, corresponding ‘to the larger 
distances between the source and receiver also provided better results. Observing the 
sound speed profiles, we identified a sound channel spanning from a depth of about 30 m 
to the bottom. With the modal decomposition, we verified that the lower modes, which 
were closer to the channel axis (~ 45 m), carried more energy. When selecting a receiver 
above this channel, we observed poor results without any match, and the best results 
corresponded to hydrophones positioned closer to the sound channel axis. Another 
important fact to note is that the depths where we obtained localizations corresponded to 
the ones containing larger amplitude modes, mainly modes two to six depending on the 
source frequency chosen to be sampled. 

Regarding the use of matched-filtered bs raw data as inputs to the FACM algorithm, 
our work suggests that it is not a major issue - one can have localization results applying 


raw data as useful as the ones from MF data. However, these results were based solely on 
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LFM chirp transmissions and may not apply to all signals generally. Further studies are 


needed to examine this issue more thoroughly. 


B. RECOMMENDATIONS FOR FUTURE WORK 

Once we have data from 30 of the 32 hydrophones in the VLA, an almost natural 
consideration would be the multiple-phone analysis of the localization algorithms. While 
computational load would certainly be an issue, the comparison with single-phone 
analysis, and possibly better results, would be worth it. In this case, one should consider 
the use of the array tilt information contained in the data files to correct the arrival angle, 
vertical structure, hydrophone depth, etc. 

Another possible approach would be the use of non-traditional processes, such as 
artificial neural networks (ANN) and minimum-variance spectrum estimations (MVSE), 
both in the arrival structure analysis and in the data preprocessing. This could surpass the 
over-simplifications of the a used techniques and lead to better results (Ma, 
Y.L., 1997). Also, since the shallow water environment is generally a dispersive channel, 
and presents effects such as strong bottom reverberation, signal distortion and noise 
fluctuation, an analysis of the signal mismatch influences could introduce improvements 


to the process. 
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Throughout this work, only three scenarios were considered. As the NUWC-TSP 
experiment provided at least 14 other sets of data where parameters such as source 
location, signal bandwidth, and/or pulse duration differed from the ones we studied, it 
could be useful to process them to try to confirm our results. Moreover, it would allow a 


comparison of the FACM performance for longer pulses and higher frequencies. 


I 





LIST OF REFERENCES 


Bendat, J. S. and Piersol, A. G., Random Data: Analysis and Measurement Procedures, 
Wiley-Interscience, New York, 1971. 


Brune, J., “Transient Localization in Shallow Water Environments,’ Master’s Thesis, 
Naval Postgraduate School, March 1998. 


Chin-Bing, S. A., King, D. B., Davis, J. A., and Evans, R. B., eds., “Lloyd’s Mirror — 
Wide Angle Propagation,” in PE Workshop II: Proceedings of the Second Parabolic 
Equation Workshop, Naval Research Laboratory NRL/BE/7181-93-0001 (US 
Government Printing Office, pp. 62-67, 1993. 


Chiu, C.-S., and Ehret, L. L., “Three-Dimensional Acoustic Normal Mode Propagation in 
the Gulf Stream,” in Oceanography and Acoustics: Prediction and Propagation Models, 
p.179, edited by A. R. Robinson and D. Lee, American Institute of Physics Press, 
Woodbury, NY, 1994. 


Chiu,C.-S., Apel, J. R., Badiey, M., Finette, S., Headrick, R., Kemp, J., Lynch, J. F., “An 
Overview of the 1995 SWARM Shallow-Water Internal Wave Acoustic Scat De 
Experiment,” JEEE Journal of Oceanic Engineering, July 1997. 


Defatta, D J., Lucas, J. G., and Hodgkiss, N. S., Digital Signal Processing: A System 
Design Approach, John Wiley and Sons, pp. 628-649, New York, NY, 1988. 


Gradshteyn, I..S., and Ryzhik, I. M., Table of Integrals, Series and Products, Academic 
Press, Inc., pp.966 and 972, 1994. 


Hager, C. A., “Modeling the performance of the Pt.Sur Hydrophone Array in Localizing 
Blue Whales,” Master’s Thesis, Naval Postgraduate School, September 1997. 


Hardin, R. H., and Tappert, F. D., “Applications of the split-step Fourier method to the 
numerical solution of nonlinear and variable coefficient wave equations,” SIAM Rev.15, 


p. 423, 1973. 


Jensen, F. B., Kuperman, W. A., Porter M. B. and Schmidt H., Computational Ocean 
Acoustics, American Institute of Physics Press, Woodbury, NY, 1994. 


de Kooter, P. M., “Variations on Autocorrelation Matching and the SIFT Localization 
Algorithm,” Master’s Thesis, Naval Postgraduate School, March 1997. 


> 


Ma, Y. L., “Emerging Issues in Underwater Acoustics Signal Processing for Shallow 
Water Environment,” in Shallow-Water Acoustics,” pp. 535-537, edited by R. Zhang and 
J. Zhou, China Ocean Press, Beijing, PRC, 1997. 


Miller, C. W., “Estimating the Acoustic Modal Arrivals Using Signals Transmitted From 
Two Sound Sources to a Vertical Line Hydrophone Array in the 1996 Shelfbreak Primer 
Experiment,” Master’s Thesis, Naval Postgraduate School, June 1998. 


Miller, J. H., Benson, J. L., Chiu, C.-S., and Smith, K .B., Transient Localization Project 
at the NPS, Technical Report, Naval Postgraduate School, April 1996. 


Pierce, D. D., “Range-Dependent Passive Source Localization using data from the 
Barents Sea Tomography experiment, Ph.D. Dissertation, “ Naval Postgraduate School, 
June 1996. 


Press, W. H., Flannery, B. P., Teukolsky, S. A., And Vetterling, W. T., Numerical 
Recipes, Cambridge University Press, New York, NY, 1988. 


Smith, K. B., Personal communication, Naval Postgraduate School, April 1996.: The 
Monterey-Miami Parabolic Equation (MMPE) model is an upgraded version of the 
University of Miami Parabolic Equation (UMPE) model (Smith and Tappert, 1993). 


Smith, K. B., and Tappert, F. D., “UMPE: The University of Miami Parabolic Equation 
Model, Version 1.0,” Marine Physical Laboratory Technical Memo 432, 1993. 


Smith, K. B., Rojas, J. G., Miller, J. H., and Potty, G., “Geoacoustic Inversions In 
Shallow Water Using Direct Methods and Genetic Algorithm Techniques,” Proceedings 
of the Pacific Ocean Remote Sensing Conference (PORSEC ’98), pp. 703-707, Qingdao, 
PRC, July 1998. 


Tappert, F. D., “The parabolic approximation method," in Lecture Notes in Physics, 
Vol.70, Wave Propagation and Underwater Acoustics, pp. 224-287, edited by J.B. Keller 
and J.S. Papadakis, Springer Verlag, New York, NY, 1977. 


Thomson, D. J., and Chapman, N. R., "A wide-angle split step algorithm for the parabolic 
equation models," JASA Suppl. 1 83, 1983. 


Tolstoy, A., Matched Field Processing For Underwater Acoustics, p. 7, World Scientific, 
River Edge, NJ, 1993. 


54 


INITIAL DISTRIBUTION LIST 


Detense Pechinicalimlormation © citer tee cee ee ae 
8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, VA 22060-6218 


DudieveKn@xLIDrary s.cs:ciseet..c.segeeeenee coe meee eee en cn ee ee ees 
Naval Postgraduate School 

411 Dyer Rd. 

Monterey, CA 93943-5101 


Profesevanl oie SING ae. ....ccsc.cen ee eee ee eee as ee anes once 
Physics Department, Code PH/Sk 

Naval Postgraduate School. 

Monterey, CA 93943-5117 


Prog, Chins -S amg: CW siccic 20s .cidecsabecstereeee cco emer oe ase tne earn eee eee ne 
Oceanography Department, Code OC/Ci 

Naval Postgraduate School. 

Monterey, CA 93943-5122 


CDR Mitchell Shipley ......2:22.:620¢5.2051 esac cce 
Physics Department, Code PH/SM 

Naval Postgraduate School. 

Monterey, CA 93943-5117 


Stephen Greineder sscssss 2. ssscc, eesrceeoe eee =2 ante en cne scenes eee otc 
Naval Undersea Warfare Center, Division Newport 

Code 2121, Bldg. 1320, Room 380 

Newport, RI 02841 


Instituto de Pesquisas de Marinliar vite eee ere 
Rua Ipiru, s/n® 

Ilha do Governador, Rio de Janeiro, RJ 

BRAZIL - CEP 21931-090 


LT Arthur F.Bettega Comea..2 secs 
IpgM - Rua Ipiru, s/n® 

Ilha do Governador, Rio de Janeiro, RJ 

BRAZIL CEP 21931-090 


as) 












12 “9° 3974 


10/99 22527-200 wre | 


1, 
ca 
z 




















in es So A 
ii 1 OUDLEY KNOL DE i Me rhe ate ot are a rea) ot 15 A 1 r 
(| NI SU UMUSV EPS DEED UD LEED + PUL yOn TL ADR ete ee ae % 
oat tl {| | } | \ i} IN| ! | | | i | iI fy 4 ¥ Me au 7 OP . ah eee 7 fa. 
ae ( MTA TALL HT | AA HI MTT | i Ae a oe Ue yee ee ca NE Pes 
, Ky | lah {I \| Mt Sat 7 =e 5 Pa 
oi int | Nt {| \| | \| | Hy | | tH) Wh it M | I ' 5 Var ; faa dN Crema eel ee f f p fi Be a 
' ae ATES EEUU EE TE iT NAVA AIHE Ace Pa LE Fev b hs Mia ee he 
a 8 00368173 5 he INO ec co 












CoC POR RL 
Py eee AeA &3 
' re | 





Ai cn 
ea 
fen a 


























Fee Ue - 

| ngye topom eee | ' ibe 

aha OI Ea Pc ee 1 ae a Pn 

ay ere ry a ie oa t us qe at o8 out Jat ' 

A > 1 nr nt ‘ 1 t i! ai ‘ 
oe F 


Pra ee 
be Te 
Nad RU ae 
a6 4% OE fet hie 
PPT ee ee fy 
av the tah 
14nd" 
















1, mats 
A ie Ay nt 
ey ere 
a 












a Ae a ms 


















ee 
Ars HA 
i af Iman rT 
fe ere Ca TR kas 
This tte ta ftrant ly eh ries Ta rae | tthe 
f Pg yin 
Te Lae 
FP 


fake et ekgaacem aba 
F’ ‘4 
¢ Cy a th 
P Or 









etre ri Pen 
WRT Nal 5 ree PRP EL et 
othe pape ee a ert i AE Ce eee a 
ze Pn: shube 
Pa: kt * ry F 
oy ; 7 i ye PP 


eer MT ark: ef a 
me Salah ag 
ets rea ped “ut P ea he “eee 
Taeweet Sr ! od rary Maa eee Pac aas 
Por Dae rea 
Te Lag 





































































































































































































































































































































































































































































































































© rere a 
Forder sath ey edith ed om] 
eRe elev RR TG f 
ye Sateen Fee OA hart ME ee rh Reet Keys A 
a PAs rt re i & z 7 ry rr oe ay ey Fi , 
ate Ree eectet Apart eis t lee na apse ab oon 
Fd fener rhe j i feria 
Pe} ge opr OA jeter oa , A, foe 
Sf roe) - OT . 4 . ha mane 
eee TN ato AA cre er LC 
eed Ci firalatatert am ne tiie Fi 
ee ee aS tphglgl 4 4if eas te c iD of Ly LU 
i Pe eT el PPE ia ee PP fe thet muses f Te ee See us ear] 
“4 Parag PES ae ee Une al cee st eee f A att be Ue Ree Aes a a 
re) Pate Ld H ft ntea dy hot Td e les } Oe OT het ay rh 2 ArT SATE aay Up OA em nates ree Ge , 
var EYL Heit vevy e'e 14 pe Fh ee Se % re pike Contec ey recie aia eto eC Sea 5 - Fl Aer, a, i ota ae i 
Poe ae 4 +. Ly Ar hana Pe oe re Cr) zn Pome 1) te 1 Pe eae eure n ry 1 ¢ f 
te a Bas feta Pat Pee a wT Par ‘Wear rae ie Po ee Pre me Pt pee at P tn Fis. ta Pa Pa 
wy Ve Prt Weherers U rt PU Ae a y t Parris le, tata oe Ser tt: Oe eae ' ieee ea 
y yl ge Cr ee Sir THe neti See ea t Baer CC Ae Ce IP ls Pr A yt " ee oro 
eee eae ea 1 ta tet aaa ! Pree ee oe) PCr ee Pree) Lee hal Aferies tke baie eee aL) fat Co ee qugeed’ Lie Mee s ' 
rer is \o4 Ae ey Lea eh eth pind Fb ert & ri ani H Seca ay rl Pa r Pi 1 Pe F A ; ry 
= ital , ' . a [ter We a ode ti e {tee se Cy 
Pirie ee ‘ ( H i eer + ne 4) ' a Tay aAeneswore i id an 0 s ‘ ff e ia) 
5 Fie SUL! te ee oe vel al bec baal f eu oad 
hed ba a ener erl eee ee rar a Ue pee » 8 
a te rips Agee Te ec Se Fy Pir PHL Pen Se ‘ 
free a Often ay Le itl é Party a PnP ee “uw er A ne 
RY Te a) tire a hd ee a) Peri feet ei) Rar Rice Wee pint Co 
Pe eeake oes ay ee Pare Y ha) Pier PC Te ec ee Pe ee A ti 
Pee) SO er es hybidpten ro rf tat wet fanny ones et DOC Ter ie icc ee ve) ar) 
PPP a cou pot ot et Pr er ee athe ort ore Ue 
TOTO a ec Aer Ta ee Pare iT eee de Te . ou 
* ot eee aa TT PPT a ee Sul ee ar eer Oe ‘ we fe m4 
PY Se CID Ck Ls a af Mere Pan LT ae ae ed Pr ae eh aes Pre aU uw 
at tr o7 Pd im CSL Peete A bet Coney aed Pee Zt ea Geel aLL Pare | Fara Oe mend eaenee er ein 
(prt 54 are Pare OTs Saat HOE UE ean dd UA Cea Te er ere ee See Peay OMNOuL) ded ir Oe a) Pr 
ae eee i) eae pis da Se ami? tat here! ae ie fe a eA t Ea Ri 3 u PUN aE ere i ee Pa Pe i 7 AT ile ete * ar) 
Fy c! eee ee Le ae 5a hai a AS Se TL een See ee in Pina ic ot | Pen aye (ea 
are ee Seis ref Pabriees! Looe POR Pie ide a“ A erg he Wd 49 & FRETS LOTTE AOeL Pik Sal Cs ca Vth ae MTLzs ROR 
; ath I 3 Ay ee Pe ee arene tm ity? ri rit srr) CataList 
mea CY a te CN an é sad A , Paar PY ae PE SE pues a hee p éafel a Scie b anata 
Pe? afaFe hie ; A a Fey , n nibsgered TH Te PrP eM i Lc a vee ‘an 
t sr ST eee a Pent ie Tet el SAE BLEU POC ie: | Pa er ar 
Pa hiner! Pre eee ea PRT 21. ae oe ie eee Bere Ti Comune dh ‘ ss » : - 
rine og EA aT ben gah at : Rees et Sy Laer ta Perey an Fi 7 - 
ries : MTP Um br Cer rs tes Pr TLL eee ee Pe Fi » 
er Ot PTT liad CU potiggesenend YT 7] ,vie : 
ee) f ‘i es “UP ra ae aif , 
n A ae ee yted nea 7) ar 
y 5: Hoe 0 a Ed s @,f ¢ mY. 
rf ‘OTA AK a) ACY Cite SU seOUi aa Kae ‘i “ee AL 4 7 te er 
Pea ee o ' CAE 3 U SO Ch eee dr el Cr rin ee) pee Oe AY ere LL 
RRO amt sl ae A ee TMS rae teed Ae pet SLC De aIDnL, Galea i. 
eres eT Le al Teh Es CL ee mmeede pe ae Peer : paeen ripe eet ute 
ie} r2ebargdard ad 4 4a a eae a TEL 1 Pi er a aCe Pree 
3 1 i aD Teas : Pr Tame Ti ye i Prd ernea to r Pa 
re Le. u's Far Lee oe Ceti erie 00° 0 
raters + ce " pe ft fe PI lait dale } an) PT Te i} oo 
rs 1S CT Cy OP te ' in P Pa 
ait te a , Ph Yoh Le Py ee > ether eo bP rageen # rn 
@yta be - rari tt, RT te ay ae ut ‘ 
1 APT Y ae eae aati 
aa -" Py Ramla d * DO 
n ' Uj 7] gr P ' « 
eee hd pye ge rate a Ted 
nore : D ry PmarThty oat) Carat tar, F 
at ref Fr " f ete Ai ae A 
"Pet ee J. ree Oe TL a Ped i) eG rl 1 
Ameer TL ee Lee ee La , e S d 1 n rh i 
a eo a Tes! pte tae lie can 5 Ye Ty ' 
aT La ad Pe ee rT ee Fy D Pie n ‘ f 
zt vc eer Dee ery Cr? fee Pres Yat Ln a ade ioe rr a} 
ia et ik Ue ote 7 =F) i er 2 i a tt @ 7 oad 
ae rae ean ot) on nh wet erent T } ? ' LU 
FAME For ARR tage 2 § co) Fl Pan) A 
ad ee a tu hg i 1s toa rd 
Sen dhes an é ry wigs hn er Fy 1 Pa F LCi he 
a re ya te a | inet tra 
coyas Pie mote ye 2 ® b wz ‘ | a ae 
Co | a es 02 fe ran PEA} a r . ae . 
oA el ss Bb oatyrst 6 Pe a i) a a | o 
"I A i Ls iy ‘ é ? 2 8 7 fn 4 rd a? Pa 
$ CB at a a aie e fanens Pe rst E |e Cee, me Fi A 
: 1 I i Pan i in ' au fy td 
: + eas er mn it ts Spee refats we re 7 1 
aN , " ’ io Sit ree ee ed A a2 et 2 i 
ny s ote J Pe char et raed aT ee ae fe i a’ ' ‘ # oat Me 
Mr ae ple ttats feel so “6 pi Ad E VU is Ai ae st te!t «# ee 
>. CA eee int ates rf es COL a ern A ay pre feet 8 Peers Yn tv 
Me Ao F) or 1 st Oe ‘ 4 ' rl ‘ yi ‘ @e of FY ‘ 
2 a) Ae yi ae 7 i ) : RT LP eae Lees r LE eo Ti carl A Jeon 
t agg ) at ¥ ¢ ad 84, \ ® b ‘ ‘ ‘ a] rf b rr J 
ed 5 i C 3 A PE hd ri re stgest F “tear F rn) ‘ 
ogee H : : UP. ar) ae rn Tr Pa) Pe F a 
Pu A ‘ Pai A one a A ‘ “a £OF 
a5 td ‘ = Pay ag *) rr ayy 4 eebor tt i FD nl F 
Ry ae | iF Ps nT) A rs | r; A fa 3? ro a 7 4 
- tld wee! by P| H , 5 Phy x, a a) iG a? i 
ng at a*5 Pho] Us rt f a 1 Pe | ad aa a L iM : 
8 te al Pie ie be Leon Pet) : "1 MU ic. o : ri rua) ‘1 . F or r 
4 | ri Tra fe F’ ate » ! Pee | ry ead e F Co) ny ‘ 
7 1 oft ws Pf x 5 a a ts oa D 
. Ss tor oo c Ty A " rr Pi a 5 ? ae " ry 
oe ‘ setae) ot ea © ef J 1. ai ks ! ra 5 ] ad F et 
° +r hs e tu te 12 @ . ry U F] Fy Py is hy id ° 1 
a , ‘ Pe « fecha ‘ oa) Pee ‘ F 18 ] { SY 0 i 
Hi ; ee ba bites Pa Loe, 5 ‘ k ‘ i » . 
cd ie ead | > vy a ry i b “hoes 5 
$ ey 7 See ee nt Herr y | r O ie ' ' r Pe 
Pit, “y Pa . Suter eee Paria ro ebe * OU 2 1) é 
A “ig A en ee iret a A we uf JP : r 
ry | ad 7 eset ai) ‘ we t « r cel 
i reed ee ee So ‘a 78 rl He t ‘I ‘ an Ld * 
x yy ° By net Set Bake Fe ac dime fe ad i r a a a ar tal OE hs ‘ 
Fg F A a eet ae sary | rs ae Pa i « , Pa ae Pn ° 
, ie Teg 1 Pam 2 V2 ee Pe -) 7 ¢ r iu - : ee y 
Wy ‘ erie be 4 o fi td Tee ee ‘ y 1 
b it pris 2 vr =f A Fy ry bir ts pt f ‘aa the 1? ' 
ae 2 tds ry , bs ry F vr 7 tN 8 8 D rie oh Ae ee Le 
Le) be J bry: a} $s ha ‘ or ee 2 Fy % ra ee 
TST5 nue ray ry A , Lf et %e# oF) , a # ° a 1 
ey 43° es J rT i , = “s r > ah rar A ‘ ) ‘ a> . ' A 
iar aly a ind "a4 ' r n Pr r ar) ' 4 ) es 
ay eee Tit Clb i ond arn cs Sk ed bord .t 1 et ct ah? ? a ‘ c f A 
{ares " : 5 ." Pr an ar Pa 4 Ul Pr RAE o n 
in tad rede. Pe Ps , Pere Us > 4 FI { f 
o LAE ] fe te ves Py n ti. ei Treks ‘ 
‘ Py ee | ae F 7) 4 i | 7 rT) F) ee Aa | ae af ¢u Pa 
Z ri fn ae ed a “yn fet Py 1 os ris t} rr O My J 
ae ees mae ot ee aD F rs Py P . , 1 
pt Pree . one " Sh e* e F ' Pr] 
te yt ete a7 y e tr F Fi n 
7 To Boe bod Peat I ig Sn y! etter yt best 8 Pt MA ‘ oi Ieee 
ad are IM) aN , i : yo? peg yae* ‘ tox $8 — 
PS a , ne A - . - 4 r 7] . € ee . A 
ae 4 CaP Y a f Pa J ‘ ty 
ihe fal r ‘peere GAS ¢ p ee pi ae “ee ’ Pa ae 
tre Pe 1 ard bs t a 
> PL Py n rel! ¥ a ck) ae oT ¥ , 
Sd si LS = .’ i 7 rh , ‘ i 
-o Pha Pe + ite eee Se A ot t ney? ‘ 00 ‘ ce Ay D VL 7 4 
Fargewerte d a PCL e 4% y beens PZ aa ° 
ee ee ae AT i oh a] AY eee Oe eof 4 H n eae Fa 
Pd ka bd rte rad rh o‘eu Pa ee os a) oy Hyiee ren, 1 ry; Fs ee . o ar] rl 
St ie eb Toe oo ty ert ea i a Ue ee | nt ee un Pe Fi * r : ; 
ies poe fds head F n PT ’ the wh nae Aa ® A ° 
Lt eS oa) » aL a eae ee Pe ee fe a oe er n f ¢ 
oe ane ry a Sahel a Ta ty, ry {sas nar ‘ ; . : ' 
© U (re 1 ' ie by \ rhe oe © tas pte ear y | a2 ‘ we ‘ Ts es 
amargernry Parl et. P +: rl ” rat ! nrc a Y Poe ee bed be : 
a on tty Me Pe er ee ied ‘i e te n Pre Sa a L & 
had s ery ge! vi Tr rae 4 ie et r 
o egeeg,see in I *{ Oman Ya Aa ey nex Th a aaa i eset ' 
Phd Peery 2 eed Pd wy Tat Pa Aa r rr Fy ‘ ors : 
cana eee a TC L 4 Le 1 ot & 7 5 ‘ U 
p c ey eed Pa nt 7 - , ‘ 
ey eeCtrie ails % } ee pi 
re rt Let Ale rt 
cher Siar p ’ fy 
Pr x s + 
ry rv ee | nr a 
af Oues oY] + ‘ tT 
ean qeerghy? ' n O ae 
fl o 4 . P 
5 ee 
y § % Chen Yon Yat — 
F 4 Fi “ 
ran : 
¥ ALT hal ek a i 
at rd hr TY 4 A ‘ 1 i 
hs rg 9 aL gu yh gos eens A i” ' 
1% Oe a! rere | year Ore atts H vi D ' 
Os as ig oF rs bs , PS be Leno) ste teahelhegee ny slewat ey teal aU ry 
LP a tebe A A A RE bd ahd | peeepen peg eats? asae tate 68 yt et ¢ a en A ‘ ‘ ‘ 
i ry Agree ye Te Pra Ret tl eee | aeovm™m 4 PPh 4) ee Py i} ira | | iM f 0 e @ ' 
Ht Fi 3p SAL aL SEA Cr eer 1 a REPS ST fie a ee Peta YL oe i A a! e 
set pret pb eer RYT aes et a Pe i ee ey ef i Parva) a papeeh toe er Pan Cn} Ld ' ‘ J 
oy rts ee eda LU sau fae g tte Mrbyryeyhe vee i es ‘ Pr Le iO r ‘ n q p 
peer aS Loa Biel marie Tepe gL Uae aI. a ar ‘ ri a | , ‘ PUL. ' , 
ee hd a | Par yy aun aie 1 bt mes ra ny gt e Uw ° ny 
sere 7 fi a Pi rt PA We ll yet oo "¢%) ” 4 ' Fs 
ry + ‘ 5 1 ape } team dee RNa ke m ni Pa r ri an] A ‘ 
wee : oT p . vy Sea elie Way A har We Seed ao ett eee var eee tee 2 
errr 7 ro hp ’ 7 Fyre PORTS Here aia ar tL tye a ie e ethan a ae 
art haba . t ) 5 ? 2 Se a ee te ay 8 =F | A . oo » e . 
at tag aE a Ay : . ari WESTIE R Ln Lie eet oe od bere fog ts LC A A a Ae ee f n 
y 4 * aon Pe ar? boll Le hh | a eee TL esis mus ae) " ‘ 
; J r . ae eC at e cette its A ee » . ch ad a 
at irk er Shido Aer ace rin ey | ere rar 
Te es Cn. marl plane ty wn ey git n i ra Fa 
(eel * tae Tb Ve ak ite Lee ee ee Patt Ty] A Rete er te ee s a * 
is Pa Leake Bal Pa a ereaay Ar Thy a 1" PP et | 4 i rl ' fu , se 
a a i aL Pa 1 eg? ns rn re a ae ' ° 
rped,t 4 ei Pr a re) peree =: Pe Ante n a) A ' ry 
apt PCM Y marie “" i 5 ry ‘ ir ar] SY e 
' te ee | ¥ 0 ‘ ‘ a! Be 
i LOC Ee. er z ; ' F re f Oo O ; 
cant LT o4e he a ine ry a ro yt i 
LI L rad oh 5 ads Bt ne ari Om he PT ae A ee h ver) Lie LL) ® es 1 : bs 
OTe Lhd ee jl HAE ene etki cate f Cacia ie ies ma D ¢ eef ‘ a 4 4 
Coe ke 1 ar ee era Xu A reget fit iL are ere y n ' 
wa) MPL a ee or eee a) Lf Pvt Uke rae Nel ati. a A v N A O Cf 
H tt Cree ery ATTY Oe li hl he ale 1 Pris ar aouwuntGet tee anf rue ge at ‘ 
Ca aaetior | oS RT gheete tpg ang erate tt Se Te c om 
Dit dle BE PRs ted etd ee ot Lie Ya er Dl gue were Pr? ‘ 5 5 e . ‘ 
OD ta v Spengiae ot rl oi th antd taht 5 Poe 8 ee Y 
Ta ; Ce eg Dey mn Pee ATL ‘ ot J bd 
aa Shand ny! Aan * Pattee ‘ A 
Tt a Ly ape | 4 F t , i 1 
OTC} hd Bee Ad | 4 Z . a. A ‘yf om La Fy a 
F 1G a 1." ' An a 
rar bal eee Ms sey sayy ele er ek Pk a | ed reat Fr +7 ni rae | rh a | rae 4 
ey : . ’ tM ue fj re Oh roa 
orpa a) oe - aE. Pt Cee recat ot : 7 
Pete ; Aye ar ee Ra het 5 x 1¢ y 
Se LU bat iol ra hr ti oo a5 WC | 
Pe ed ‘ 1" Pas veer ray “a ire s ¢ . 
so Tera) oo ae er ey ae ILD Ten Ca akat a : is 
A ey ee deren! Ae aoe 4 
eee LL Td a bk ‘ 
ran ee & ee etate r ‘ 
pet phe u A a) a : ae 
. ed ae | u Pn Pn rd) 
A A eee Ce oe i 2 
poems «*4" 1 it tary Pa a? a 
Li . Ce a Fo ene eke A y 7 rs ‘ ry 
rs TS. r ony oeme ot f be ai 
Tiers ee of phaad ve ae ee Pa 4 i 
ee eran re ee oa o : 
aia ott athyerts bt ree.) 48 
MCS TTE th ae Ro ' 5 1 i r 
~~ nt A a 
1 gf fa ' a 
ft 1 Ve hy 
x 
n 
Pn 
1 


